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ABSTRACT 

The analysis of the multicolour photometric observations of MW Lyr, a large modula- 
tion amplitude Blazhko variable, shows for the first time how the mean global physical 
parameters vary during the Blazhko cycle. About 1 — 2% changes in the mean radius, 
luminosity and surface efii'ective temperature are detected. The mean radius and tem- 
perature changes are in good accordance with pulsation model results, which show 
that these parameters do indeed vary within this order of magnitude if the amplitude 
of the pulsation changes significantly. We interpret the phase modulation of the pul- 
sation to be a consequence of period changes. Its magnitude corresponds exactly what 
one expects from the detected changes of the mean radius assuming that the pulsation 
constant remains the same during the modulation. Our results indicate that during 
the modulation the pulsation remains purely radial, and the underlying mechanism is 
most probably a periodic perturbation of the stellar luminosity with the modulation 
period. 

Key words: stars: horizontal branch - stars: variables: other - stars: individual: 
MW Lyr - stars: oscillations (including pulsations) - methods: data analysis - tech- 
niques: photometric 



1 INTRODUCTION 

Neither multicolour photometric nor spectroscopic observa- 
tions have been previously obtained with good enough time 
coverage of both the pulsation and modulation cycles of a 
Blazhko RR Lyrae star that they could be used to derive 
any phenomenological conclusion about the changes in the 
global physical properties of the star during the modulation 
cycle. To cover the pulsation variations of different shape 
in each phase of the Blazhko modulation with observations 
several hundreds of hours of observing time were needed. 
This can be achieved only with telescopes 'dedicated' to the 
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study of the phenomenon. As 80 — 90% of the telescope time 
of the automated 60cm telescope of the Konkoly Observa- 
tory is allocated to study RR Lyrae stars, during recent 
years we could first obtain BVRcIc photometric time se- 
ries of Blazhko variables that are extended enough for such 
an investigation. 

We have detected 0.010 — 0.005 mag systematic changes 
in the intensity weighted mean (V), {B) — {V), and {V) — {Ic) 
brightness and colours of RR Gem and SS Cnc in differ- 
ent phases of the Blazhk o mo dulation (see Fig. 12 and Fig. 
9 in lJurcsik et al.l (|2005l ) and lJurcsik et al.l (|2006t ). respec- 
tively). The intensity weighted quantities indicated slight 
brightness and temperature increases in RR Gem at the 
time of the largest amplitude of the pulsation, while in 
SS Cnc the brightest mean magnitude and bluest mean 
(V) — (Ic) colour occurred during the decreasing pulsa- 
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tion amplitude phase of the modulation. The small ampli- 
tude of the modulation of these stars, and also the ani- 
bifiuity of the equivalent stat i c colo urs of RR Lyrae stars 
l|Bono. Caputo fc Stellingwerj Il995l ) make, however, these 
results somewhat uncertain. 

The light curve analysis, utilizing mostly the V data of 
our ~ 1000 hours extended multicolour CCD observations of 
MW Lyr, an RR Lyrae star sh owing large amplitud e Blazhko 
modulation, was presented in I Jurcsik et al.l (|2008l . hereafter 
Paper I) . The large amplitude of the modulation of MW Lyr 
makes it possible to detect changes in the mean magnitudes 
and colours more definitely if there are any indeed. 

In order to exploit the most information possible from 
the multicolour light curves of RR Lyrae stars we have 
recently developed the IP method ( Inverse Photometric 
method. ISodor. Jurcsik fc Szeidlll2008l ). This method gives 
good estimates of the mean physical parameters and their 
variations with pulsation phase exclusively from photomet- 
ric data without any spectroscopic observations. Using the 
extended BVIc time series of MW Lyr the IP method shows 
the difi'erences between the pulsation in different phases of 
the modulation, and whether or not there are any changes 
in the mean global physical parameters of the star during 
the Blazhko cycle. 



2 DATA AND METHOD 

The photometric data utilized in this paper were pub- 
lished in Paper I. The observed colours and magnitudes are 
dereddened using E( B — V) = 0.1 mag interstellar ex- 
tinction given by the ISchlegel et al.l (|l998l ) catalogue and 
Av = 3.1E{B - V) and E{V - Ic) = 1.27E{B - V) rela- 
tions. 

In the present paper the mean light and colour curves 
of MW Lyr and the actual light and colour curves in 20 
different phase bins of the 16.546 d modulation cycle are 
investigated. The colour-magnitude and two-colour loops of 
the pulsation in different phases of the modulation are shown 
in Fig. [1] This is the first multicolour photometric observa- 
tion of a Blazhko variable that makes a detailed study of the 
colour behaviour changes during the Blazhko cycle possible. 

The analysis is perform ed using the IP method 
jSodor. Jurcsik fc Szeidi 120081 ). that finds which T^// and 
Vrad variations during the pulsation results in luminosity, ra- 
dius and effective gravity changes that best match the light 
curves using synthetic magnitude s and bolometric correc - 
tions of static atmosphere models (|Castelli fc Kuruc3l2003l ). 
The input data of the IP method are the Fourier fits of the 
time series data. Solutions are obtained using the Fourier 
representations of either the V , B — V , and V — Ic ot the 
V, B, and Ic time series. 

Those parameters that definitely must not change with 
the Blazhko period, namely the metallicity, mass and the 
distance of the star are determined from the mean light 
curves of the entire data set and are fixed when studying 
the changes during the Blazhko cycle. Two mean data sets 
are investigated to determine these parameters, one corre- 
sponding to the original observations, and another where 
data are corrected for the phase modulation with appropri- 
ate time transformation of the observations. In Paper I it 
was demonstrated that if phase modulation is due to pulsa- 



tion period changes during the Blazhko cycle, this treatment 
is valid and the mean light curve of the time transformed 
data represents the mean pulsation variation better than the 
mean light curve of the original data. 

As no spectroscopic observation of MW Lyrae has 
ever been obtained, the metallicity is calculated from the 
ljurcsik fc Kovac3 (|l996l ) formula which derives the [Fe/H] 
from the period and the $31 phase difi'erence of the V light 
curve. This formula gives —0.6 and —0.4 [Fe/H] values for 
the mean light curves of the original and the time trans- 
formed data, respectively. Therefore, it seems to be a reason- 
able choice to use the IP method with [M/H]=— 0.5 model 
atmosphere grid for MW Lyr assuming solar scaled chemical 
composition distribution as significant enhancement of the 
a elements are observed ty pically only in more me tal poor 
variables (see e.g.. Fig 9. in IClementini et al.lll995[ ). 

The mass and distance of MW Lyr are derived by apply- 
ing the IP method using the [Fe/H]=— 0.5 model atmosphere 
grid on the mean light and colour curves of the entire origi- 
nal and the time transformed data sets. One basic input of 
the IP method is an initial Vradif) curve (75 means the pul- 
sation phase) which is v aried in the minimization process. 
ISodor. Jurcsik fc Szei H (|2008l ) defined this initial Vradi^p) 
curve twofold, first using the template Vradif) curve scale d 
according to the Ay^^^ — VAmp relation given bv lLiul (|l99lf ). 
secondly utilizing the Vradi'fi) — I c('P) relation that was 
shown to be valid for RRab stars in lSodor. Jurcsik fc Szeidi 
(|2008l ). For variables showing no or only small amplitude 
light curve modulation the two types of initial Vrad{^) func- 
tions led to similar final results. The Ay^^^ — VAmp rela- 
tion of individual Blazhko variables during the modulation 
cycle is, however, steeper than that was defined by iLiul 
|l993), and from the K-ad(v5) — Ici'p) r elation valid for non - 
modulated RRab stars (see Fig. 2 in I Jurcsik et al]|2002l ). 
Blazhko variables obey the same amplitude relation as nor- 
mal RRab stars at around their maximum amplitude phase 
(| Jurcsik et al.|[200^ . As a consequence, the amplitude of the 
mean radial velocity variation of a Blazhko star is supposed 
to be smaller than that of a nonmodulated star with photo- 
metric amplitude similar to the amplitude of its mean light 
curve. The slope of the AAy^^^ / AA{V) amplitude ratio is 
30 for regular RRab stars while it is about 41 for Blazhko 
variables. The amplitudes of the initial radial velocity curves 
used in the analysis of MW Lyr, calculated either from the 
Ic light curves or Liu's template, are therefore scaled differ- 
ently from normal RRab stars. Both mean light curve sets 
are modelled assuming three values for the A Ay^^^/ A VAmp 
ratio, 37, 41, and 45, respectively, and with zero points of the 
"^Vrad ~ ^Amp relation determined from matching the maxi- 
mum amplitude value of the Vrad curve to the Av^^^ — VAmp 
relation of nonmdulated RR Lyrae stars. 

As there is no information on the changes of the actual 
shapes of the Vrad{^p) curves during the Blazhko cycle we 
do not know whether any of the two approximations of the 
radial velocity variation that the IP method uses is valid 
for Blazhko variables. In the first method when using Liu's 
template, the shape of the initial Vrad curve is always the 
same only its amplitude is scaled, while in the second one 
when calculating the Vradi^p) curve from the actual Ic light 
curve, it varies in accordance with the variations of the light 
curve's shape. Having no a priori knowledge on the Vradi'fi) 
of a Blazhko star, in order to make even less restrictions on 




Figure 1. Vvs.B — V, Vvs.V — Iq, and B — Vvs.V — Ic data of MW Lyr are shown in 20 different phase bins of the Blazhko modulation. 
The drawn curves connect the mean values of the magnitudes in 20 phase bins of the pulsation cycle. 
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Table 1. Global mean physical parameters of MW Lyr 















d [pc] Mv L/Lq 
uncertainties 




R/Rq 


log Qstat 




time transformed data 




= 37 










3595 0.74 
147 0.09 


44.7 

3.6 


6892.5 
2.5 


4.64 
0.19 


2.982 
0.010 


0.76 
0.08 


AAy^^,/AA(V) 


= 41 










3460 0.83 
145 0.09 


41.4 
3.4 


6891.9 
2.5 


4.47 
0.18 


2.973 
0.010 


0.69 
0.07 




= 45 










3343 0.90 
138 0.09 


38.6 
3.2 


6890.3 
2.4 


4.32 
0.18 


2.965 
0.010 


0.63 
0.06 


original data 


Aylv,„,/AA(y) 


= 37 










3367 0.89 
106 0.07 


39.0 
2.4 


6889.4 
2.7 


4.35 
0.14 


2.966 
0.009 


0.64 
0.05 


AAy^^,/AA(y) 


= 41 










3223 0.98 
111 0.08 


35.7 
2.4 


6888.0 
2.3 


4.16 
0.14 


2.956 
0.008 


0.57 
0.05 




= 45 










3104 1.06 
109 0.08 


33.2 
2.3 


6887.0 
2.3 


4.01 
0.14 


2.947 
0.009 


0.52 
0.05 



its shape, in the IP process an even smaller weight of the 
initial \4ad(<p) curve is give n than that has been u s ed for 
unmodulated RRab stars in ISodor. Jurcsik fc Szeidll (|2008l ) 
to allow a larger variance of its Fourier parameters. 

Table [T] summarizes the mean global parameters of MW 
Lyr using the original and the time transformed data and 



3 possible values for the A\ 



VAmp relation of Blazhko 



variables. The estimated uncertainties are the rms scatter of 
the results running the IP method using 8 settings: 

(i) Initial Vradif) calculated from Liu's Vrad template 
curve or from the Vradi'p) — Ic[v) relation. 

(ii) Large and small weights are given to the initial 
Vrad(v) function, the first keeps it close to its initial shape, 
the second allows the amplitudes and phases of the lower 
Fourier c omponents of the Vrad(^) cur ve to vary substan- 
tially. In ISodor. Jurcsik fc Szeidi (|2008t l it was shown that 
giving small weights to the initial Vrad{v) curve leads to 
unreliable shape of the solution Vradiffi) curve. The results 
of the IP method are the most reliable for unmodulated 
RRab stars if large and medium weights of the initial K-ad(^) 
curve are given. However, as we do not know how the ra- 
dial velocity curve of a Blazhko variable vary, we decided 
to use small weights of the initial Vradiv) curves instead of 
medium weights. Any result that is independent from this 
choice is practically independent of the exact shape of the 
radial velocity variations, and as so, is a robust solution for 
the detected changes during the Blazhko cycle. 

(iii) Input data are the Fourier fits of the V , B ~ V , and 



— 7c or the V, B, and 7c time series. These input data 
are not identical representations of the observations because 
the differences between the Fourier fits of the magnitudes do 
not reproduce exactly the Fourier fits of the colour indices. 

As mentioned earlier and discussed in Paper I, we re- 
gard the time transformed data that corrects the phases of 
the pulsation to eliminate the phase modulation/pulsation 
period changes to be a more reliable representation of the 
data than the original time series. Therefore, in the detailed 
analysis of the light curves in different phases of the modu- 
lation, the distance and mass of MW Lyrae are fixed to the 
mean values of the results obtained for the time transformed 
data assuming AAy^^^ / AVAmp ~ 41, namely to 3460 pc 
and 0.69 DDT© values, respectively. This mass is somewhat 
larger than the mass of higher metallicity horizontal branch 
stars in the instability strip ac cording to stellar ev olution- 
ary models (e.g., ACS Survey, Potter et al] l2007l : Padova 
evolutiona ry database, ISalasnich et al.ll2000l : Y2 evolution- 
ary tracks. lDemaraue et al.ll2004h . The IP method calculates 
the mass of the star from the pulsation equation of funda- 
mental mode variables. The masses of some of the unmod- 
ulated RRab stars derived by the IP method are too large 
on evolutionary grounds, but the combination of the results 
of direct Baade-Wesselink analysis with the pulsation mass 
formula may also lead to too large mass values compared to 
the evolutionary values in s ome cases. E.g., 0.87 and 0.85 
dJtQ masses were derived in ISodor. Jurcsik fc Szeidll (|2008l ) 
for WY Ant and UU Get, respectively, from the pulsation 
equatio n and the results of the Baade-Wesselink analysis of 
iKovacsl (|2003l ). To resolve this discrepancy is, however, be- 
yond the scope of the present paper. 

It is important to note, however, that when running 
the IP code using any of the physical parameter combi- 
nations listed in Table 1, similar solutions for the varia- 
tions of the global physical parameters with the Blazhko 
phase are obtained. It is found that only the mean values 
of these parameters depend on which mean light curve so- 
lution is accepted. All the conclusions of the next sections 
remain unchanged if other values of the fixed parameters, 
[Fe/H], distance and mass, were selected, and analyses are 
performed using the original data with other possible values 
of the AAv^^^ / AVatup amplitude ratio. 



3 CHANGES IN THE PHYSICAL 

PARAMETERS DURING THE BLAZHKO 
CYCLE 

The detailed analysis of the light curves in different phases 
of the modulation is carried out by running the IP code 
using atmosphere models with [Fe/H] =— 0.5 on 20 segments 
of the time transformed photometric data with the 8 settings 
specified in (i), (ii), and (iii) points. The distance and mass 
of MW Lyr are fixed to the 3460 pc and 0.69 971© values, 
respectively. 

The IP method fits the light and colour curves by find- 
ing the appropriate changes in the surface radius and tem- 
perature of the star. The fitting accuracy is about 0.02 mag. 
Just for comparison, in Paper I we have shown that due 
to some stochastic/chaotic behaviour of the modulation the 
best Fourier fit to the complete V data set reproduces the 
observations with only 0.02 mag accuracy. The Fourier fits of 
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Figure 2. Residual V , B — V , and V — Ic data (crosses) are 
plotted for tiie largest and smallest amplitude phases of the mod- 
ulation, for Blazhko phases 20 and 13, respectively. The mean 
pulsation curves in both phases are subtracted. The lines show 
the deviations of 8 variant fits of the IP method from the mean 
pulsation curves. 



the light curves of the time transformed data in the 20 phase 
bins of the modulation give, on average, also 0.02 mag rms 
residual, which is significantly larger than the observational 
inaccuracy. This is partially due to modulations detected 
in Paper I with periods that are non-integer multiplets of 
the main modulation period, and to small but nonperiodic 
perturbations of the modulation. Fig. [2] plots the residual 
V,B — V, and V — Ic data in the largest and smallest am- 
plitude phases of the modulation, in Blazhko phase 20 and 
13, respectively. The deviations of the fitted curves from the 
mean light curves of the data are also drawn. The fits seem 
to match the light curves within the required accuracy, with 
the largest deviations at around the rising branch of the light 
curves (phase 0.85-1.00) where the observed colours of the 
dynamic atmosphere of RR Lyrae stars definitely cannot 
be fitted accurately with static atmosphere model results. 
Moreover, in this phase of the pulsation the modulation is 
not strictly regular (as it was shown in paper I), which may 
also explain the larger residuals of the solutions of the IP 
method here. 

In Figs. 1 3171 the results of the IP method for the 
Teff,L/LQ,R/RQ,Vrad, and logge// changes are shown in 
two phases of the Blazhko modulation, one at the largest 
and another at the smallest pulsation amplitude. The upper 
panels show when the initial Vrad curve is calculated ac- 
cording to the Vrad — Ic relation while bottom pane l s show 
results using the template Vrad curve defined bvlLiul (|l99lf ). 
Solid lines show the results when the Vrad curve is kept close 
to its initial shape, while the dashed lines indicate how the 
results change if the lower Fourier parameters of the initial 
Vrad curves are allowed to vary. Each setting is applied using 



either the fits to the B, V, and Ic or the V, B — Vand V — Ic 
data, however, these differences have only very marginal ef- 
fect on the results. 

Fig. |S]displays the light curve variation of MW Lyr in 20 
phase bins of the modulation according to the time trans- 
formed data. For comparison, in Figs. I9l ll3l the variations 
of the physical parameters derived form the IP method are 
shown. The left and right panels show if the initial Vrad 
curves are calculated from the Ic light curves or from Liu's 
template, while top and bottom panels are for large and 
small weights of the initial Vrad curves, respectively. 

There are only minor changes in the luminosity and sur- 
face effective temperature curves depending on the choice 
and weights of the initial Vrad curves. The amplitude of the 
temperature and luminosity variations during the pulsation 
of the star are significantly different in the different phases of 
the modulation. The temperature and luminosity at maxi- 
mum amplitude vary between 6300 and 8800 K and between 
26 and 95 Lq, while at minimum amplitude only between 
6500 and 7700 K, and between 31 and 58 Lq, respectively. 
The changes in the amplitude of the temperature and lu- 
minosity variations are about 50%! These results are very 
robust, the ranges of the detected changes in the variations 
of these parameters in different phases of the modulation 
are independent of the choice and alterations of the initial 
radial velocity curve. 

Substantial changes in the radius, radial velocity and 
surface effective gravity curves can be also seen. It is im- 
portant to note that in those solutions when the shape of 
the initial radial velocity curve is allowed to vary substan- 
tially the results are very similar if Liu's template is used or 
the initial Vrad curve is calculated from the Ic light curve. 
The radial velocity and, as a consequence, the radius varia- 
tions are somewhat different if the Vrad curve is kept close 
to its initial template, or it is allowed to vary significantly. 
In the first case significant changes in the amplitude of the 
radius variation during the Blazhko cycle can be detected, 
while in the second case not the amplitude but the shape of 
the radius curves vary. The only radial velocity observations 
of a Blazhko variable with good enough phase coverage of 
both the pulsation an d the modulation cycles we re recently 
obtained for RR Lyr (IChadid fc Chapellied[200^ 'l. Fig. 9 in 
IChadid fc ChapeUierl (|2006l ) shows that the radial velocity 
curves of RR Lyr remain smooth during the modulation and 
there is no sign of a double wave shape similar to what our 
results indicate if the shapes of the radial velocity curves 
are allowed to differ significantly from their initial templates 
during the large amplitude phase of the modulation. Most 
probably the real variations of the radial velocity curves are 
close to those shown in the upper panels of Fig. [12] and not 
those in the bottom panels, i.e, the radial velocity curves 
of Blazhko variables differ from the Vrad templates of un- 
modulated RRab stars mostly in their amplitudes and not 
in their shapes. If this is indeed the case, then there is also 
about 50% change in the amplitude of the radius variation 
of the pulsation during the Blazhko cycle. 
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pulsation phase pulsation phase 

Figure 3. T^ff variations in the large amplitude (left panels) and 
small amplitude (right panels) phases of the modulation. Top and 
bottom panels show results if the initial Vj-ad curve is defined from 
the Vr ad — Ic relation and from the template curve given by iLiul 
lll99lh . respectively. Four curves are plotted in each panel, giving 
large (solid lines) and small weights (dashed lines) to the initial 
Vj-ad curves and using the V, B — V , and V — Ic or the S, V, 
and Ic curves as input data. These latter curves, however, highly 
overlap in each of the plots shown in Figs. 131 71 



Blazhko phase 20 Blazhko phase 1 3 




0.5 0.75 1 1.25 0.5 0.75 1 1.25 



pulsation phase pulsation phase 

Figure 4. The same as Fig. |3] but for the luminosity variations. 

4 CHANGES IN THE MEAN PHYSICAL 
PARAMETERS DURING THE BLAZHKO 
CYCLE 

From photometric data different types of mean magnitudes 
and colours can be derived: average values of the intensities 
or the magnitudes or mean colours calculated from the dif- 
ference of the mean magnitudes or from colour data. The 
intensity averaged magnitudes are brighter than the magni- 
tude averaged ones, while for the colours the (mi — -012) > 
{mi — 1712} > {mi) — (m2) relation holds. It w as shown, 
however, bv lBono. Caputo fc Stellingwerj l|l995h that none 
of these mean magnitudes and colours match the equivalent 
static values of RR Lyrae stars in the full ranges of their 
possible parameter domain. Therefore, to draw any conclu- 
sion for the changes of the mean physical parameters during 
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Figure 5. The same as Fig.[3]but for the radius variations. 
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Figure 6. The same as Fig. |3] but for the radial velocity curves. 
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Figure 7. The same as Fig.[3]but for the effective surface gravity 
variations. 
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Figure 8. The time transformed V liglit curves of MW Lyrae in 
20 phase bins of the modulation. 
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Figure 9. Luminosity variations in 20 phase bins of the modula- 
tion. Left and right panels show results using the transformation 
of the Iq light curves and Liu's template as initial radial velocity 
curves, respectively. In the top panels the Vrad curves are kept 
close to their initial curve, while in the bottom panels the initial 
Vrad curves ar allowed to vary significantly in the fitting process. 
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Figure 10. The same as in Fig|9]but for the temperature varia- 
tions in 20 phase bins of the modulation. 
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Figure 11. The same as in Fig|9]but for the radius variations in 
20 phase bins of the modulation. 
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Figure 12. The same as in Fig|9]but for the radial velocity curves 
in 20 phase bins of the modulation. 
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Figure 13. The same as in Fig [9] but for the effective surface 
gravity variations in 20 phase bins of the modulation. 
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the Blazhko cycle from the changes of the mean colours and 
magnitudes is somewhat ambiguous. 

In the left panels in Fig. [14] the observed mean data 
derived from the photometry of MW Lyr are shown for 20 
Blazhko phase bins. The period changes correspond to the 
pulsation period variation determined in Paper I from the 
phase of the /o pulsation frequency. These plots show def- 
initely that each of the mean colours and magnitudes vary 
with 0.01-0.02 mag amplitude, but these changes are of the 
opposite sign for the different averages, making any firm 
conclusion of the mean temperature and luminosity changes 
impossible. 

However, after finding the appropriate pulsation curves 
of the physical parameters in different phases of the mod- 
ulation with the IP method we can derive their mean val- 
ues directly. In the middle and right panels in Fig. [TJ] the 
mean physical parameter changes are shown if the Vrad — Ic 
relation and Liu's template are used to define the initial 
Vrad curves, respectively. For each phase bin four data of 
the mean parameters are plotted corresponding to the re- 
sults with different settings of the IP method as shown in 
the panels in Figs. 1 3171 

Fig. [14] shows unambiguously that during the Blazhko 
cycle the mean physical parameters do indeed vary. The am- 
plitudes and signs of these variations are independent of the 
choice of the initial Vrad curve and of which setting of the 
IP code is used. We have also tested whether there is any 
change in these variations if any other possible combinations 
of the mass and distance of MW Lyr as listed in Table [1] are 
adopted, or if the original data are used instead of the time 
transformed ones, or if different AAv^^^ / ^AiV) ratio is ap- 
plied. Any of these possibilities have only minor influence 
on the results shown in Fig. 1141 only the parameters' ranges 
are shifted corresponding to the data given in Table [1] 

We can therefore sate that, for the flrst time, we suc- 
ceeded in reliably detecting changes in the mean global phys- 
ical parameters of a Blazhko variable during the Blazhko 
cycle. 1 — 2% changes in each of the parameters are evi- 
dent. MW Lyr is larger by about 0.047?q, more luminous by 
about I.OOL0 and cooler by about 50 K in the large pul- 
sation amplitude phases of the modulation as compared to 
its small pulsation amplitude phases. Though we cannot ex- 
clude the possibility that the detected changes in the mean 
global parameters arise from a complicated conjunction of 
different unknown processes, the most straightforward con- 
clusion is that they simply reflect real changes in the mean 
global physical parameters. 

The relative radius variation AR/R is 0.009 which leads 
to AP/P — 0.014 relative period change as a consequence of 
the AP/P « 3/2 X AR/R pulsation relation. The detected 
changes in the pulsation period is AP/P — 0.016. This good 
agreement between the period changes measured directly 
from the phase modulation and calculated from the radius 
variation proves that both the measured radius changes and 
our treatment to regard the phase modulation of the pulsa- 
tion as pulsation period changes should be real. The phase 
relation between the radius and period changes are not syn- 
chronized, the period is about the longest when the mean 
stellar radius a nd the pulsation amplitude is the smallest. 
IStothersI (|2006l ) found a similar phase connection between 
the pulsation amplitude and pulsation period changes in RR 
Lyrae. 



Our results also reveal that the intensity mean V mag- 
nitude reflects the luminosity changes and the {B — V) and 
{B — V), or {V — Ic) and {V — Ic) mean colours repro- 
duce the temperature changes correctly, notwithstanding the 
large amplitude variation of the pulsation. Consequently, 
these values can be regarded as good representatives of the 
equivalent static values of RR Lyrae stars. The tempera- 
tures derived from (B) — (V) or {V} — (Ic) colours may 
differ from these temperatures by 100-250 K, the difference 
is the largest when the amplitude of the pulsation is large. 



5 COMPARISON WITH PULSATION MODEL 
RESULTS 

Assuming that the variation of MW Lyr can be fully de- 
scribed by radial pulsation, some hints on the character of 
the variation of the mean physical parameters can be drawn 
from numerical pulsation calculations. Since there are no 
nonlinear pulsation models that can reproduce the Blazhko 
phenomenon, such calculations can reveal only partly the re- 
lations found in this paper. For our calculations we used the 
Florida-Budape st pulsation code w ith standard RR Lyrae 
parameters (see lKollath et al]|2002l ). Even the state of the 
art pulsation codes contain only a relatively simple model of 
convection. In these models only one dynamical equation is 
used for the turbulent energy, with additional source func- 
tions of turbulent energy and convective flux. According to 
our experience, the interaction of this simple turbulent equa- 
tion with the hydrodynamics and radiation transfer does 
not give rise to any time-scale or modulation that can be 
connected to Blazhko cycles. If turbulence plays a n impor- 
tant r ole in Blazhko phenomenon, as suggested bv lStotherj 
(|2006l ). then to model the effect, additional physics, e.g. with 
a more complicated turbulent-convection formalism, or with 
the inclusion of interaction with magnetic field, should be in- 
cluded in the codes. In the lack of any more sophisticated 
model only some aspects of the nonlinear dynamics of pul- 
sation can be tested in comparison with the observations. 

The natural timescales that come from numerical cal- 
culations are the pulsation period and the e-fold time of 
amplitude variation during the onset of pulsation (the last 
one is related to the linear growth rate of pulsation). Since 
the growth rate normalized with the pulsation frequency, is 
of the order of a few percent, the e-fold time of pulsation 
growth or decay is of the order of the Blazhko period. It 
means that any change in the stellar structure can be re- 
alized in pulsational changes on the period of the Blazhko 
cycle. 

As a first test, nonlinear pulsation calculations are ini- 
tiated with a small velocity perturbation with the corre- 
sponding eigenvector of the linear stability analysis. Then 
the onset of pulsation occurs through a transient, with the 
growth of pulsation amplitude. We used this phase of the nu- 
merical calculations to estimate the dependence of the mean 
physical quantities on the amplitude. 

A more physical test is to check the effect of the varia- 
tion of turbulent parameters on the mean physical parame- 
ters. From the 6 turbulent (a) parameters of the Florida- 
Budapest code we selected the efficiency of eddy viscos- 
ity (the dissipation due to turbulence), since this param- 
eter (a^) does not affect the static structure and the linear 
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Figure 14. Observed mean parameters (Amp(V), pulsation period, intensity and magnitude averaged V brightnesses, and three types 
of average values of the B — V and V — Ic colours) are shown in the left panels. Middle and right panels show the variations of calculated 
mean parameters if initial Vrad curves are defined according to the Vrad ~ Ic relation, and Liu's template V^ad curve, respectively. For 
each phase bins the values of 4 solutions with different settings are shown to illustrate the inherent uncertainty the the results. 



growth-rates of the pulsation modes. With this test the am- 
plitude of the relaxed pulsation can be changed and it again 
gives the dependence of mean values of some of the global 
parameters on the amplitude. From a finite amplitude (limit 
cycle) model the hydrocode can be restarted with a slightly 
modified eddy viscosity parameter, which gives an efficient 
way to get the pulsation model for different values of a^, 
and finally the averaged physical quantities as the function 
of amplitude. 

The two tests gave basically identical results, prov- 
ing that the amplitude dependence of some of the aver- 
aged quantities is not a transient event, but a real dy- 
namic effect. We calculated several models of both fun- 
damental and overtone pulsations. As expected, the mean 
radius of the star increases with the pulsation amplitude 
according to the results of both tests. We found that ap- 
proximately AR/Ro ^ 0.015Ayl_Boi holds. The mean ef- 
fective temperature decreases with increasing amplitude: 
ATeff ~ 100 ■ AAboI- These results agree well with the ob- 
served dependencies as shown in Fig llSI The figure compares 
the observed relative changes in mean temperature and ra- 
dius as a function of the observed bolometric amplitude (cal- 
culated from the luminosity curve solution of the IP method) 




-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 

delta Ab„| delta Abo\ 



Figure 15. The observed changes of the mean temperature (left 
panel) and relative mean radius (right panel) variations as a func- 
tion of the bolometric amplitude variation of MW Lyr in differ- 
ent phases of the modulation (dots) are compared with pulsation 
model predictions (lines). 



with model predictions. The measured amplitude dependen- 
cies are ATe// 80 • AAboi and AR/Ro ^ 0.012A^i3oi. 

We have to note, however, that the calculations reveal 
only one of the mechanisms that should be taken into consid- 
eration to understand the Blazhko effect, namely the stan- 
dard dynamics of radial pulsation. For example, if one con- 
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siders the standard relation between AP/P and AR/Ro, the 
order of the observed variations of these quantities agrees 
with each other as shown in Sect. 4. But since the whole 
mechanism is more complicated, it is not astonishing that 
the period and radius variations are not in phase as expected 
from the simple relation. Since in these tests the static struc- 
ture of the star has been fixed, it is not expected that all the 
observed variations during the Blazhko cycle could be esti- 
mated. First of all, this simplified dynamic system has only 
a very limited effect on the pulsation period (of the order 
of 10~* ). It indicates that structural changes are necessary 
to reproduce all the observational data, e.g. by varying the 
other turbulent parameters in the code. 

Another deficiency of our simple model is that since 
there is no built in mechanism besides the kinetic energy of 
pulsation in the code, which results in time varying storage 
of radiative energy on a longer time scale than the pulsation, 
the mean luminosity changes cannot be compared to the 
observations. Moreover, dynamical effects can also play an 
important role in the complete picture of the phenomenon. 
This effect, however, cannot be investigated without a self- 
consistent model which naturally provides the periodic or 
quasi periodic modulation of turbulent structure. To get a 
consistent model of Blazhko phenomenon, based on the dy- 
namics of turbulent convection, more sophisticated models 
and lot more computational efforts are needed. The goal of 
this paper cannot be to find the solution of the long standing 
problem of Blazhko effect. However, it turns out, that some 
of the observational results can be well estimated even with 
our simple treatment. Whether it is a fortunate coincidence 
or the indication that variations on the turbulent convective 
structure of RR Lyrae stars are an important ingredient of 
Blazhko phenomenon, should be answered in further studies. 



6 COLOUR DEPENDENT 

CHARACTERISTICS OF THE 
MODULATION AND PULSATION 
FREQUENCIES 

The detected changes in the mean global physical param- 
eters during the Blazhko cycle prove that the physically 
meaningful frequency of the modulation is fm , the frequency 
of the modulation itself, and not any of the side lobe frequen- 
cies appearing close to the pulsation frequency or its har- 
monics in the Fourier spectrum of the light curve. A similar 
conclusion was already drawn in I Jurcsik et al.l (|2006h based 
on the colour behaviour of the amplitudes and the phases 
of the pulsation, modulation side lobe, and the modulation 
frequencies of RR Gem and SS Cnc. It was found that while 
the modulation side lobe frequencies have similar amplitude 
ratios and phase differences in the different colours as the 
pulsation frequencies have, the modulation frequency has 
discrepant values. 

Table [2] summarizes the amplitude ratios of the pulsa- 
tion and modulation frequencies of MW Lyr calculated from 
the Fourier parameters given in Table 5 in Paper I. Only 
frequency components with V amplitudes larger than the V 
amplitude of the fm frequency are considered. For MW Lyr 
the phase differences of the fm frequency component do not 
show discrepant behaviour but its A{B) /A{V) amplitude ra- 
tio is about 0.1 larger, and its A{V) /A{Ic) amplitude ratio 



Table 2. Comparison of the colour behaviour of the pulsation 
and modulation side componets with the properties of fm 



frequencies 


A{B)/A{V) range A{V)/A{Ic) range 


fm 


1.379 


1.400 




1.263 1.185/1.336 


1.522 1.438/1.625 


kfo ± /m 


1.302 1.251/1.352 


1.577 1.505/1.623 


average value for 6 components 




average value for 10 components 





is 0.1 smaller than the amplitude ratios of the pulsation and 
modulation side lobe frequency components. 

Based on, on the one hand, the similarity between the 
colour dependence of the amplitudes of the modulation side 
lobe frequencies and that of the pulsation frequencies and, 
on the other hand, on the discrepant amplitudes of the mod- 
ulation frequency itself, we conclude that the side lobe fre- 
quencies are most probably combination frequencies of the 
pulsation frequencies (kfo) and the modulation frequency 
(fm)- As being combination frequencies they inherit the 
properties of their larger amplitude component, namely that 
of a pulsation frequency component. In this case, the inde- 
pendent frequency of the modulation is, in fact, the modu- 
lation frequency, fm- 



7 CONCLUSIONS 

As a conclusion, based on the results shown in Paper I and 
in this paper we can interpret the Blazhko modulation of 
MW Lyr as follows. 

• The frequency that characterizes the Blazhko modula- 
tion of RRab stars is the modulation frequency, fm, as we 
have found that the observed modulations are governed by 
global physical changes that also occur with the Blazhko pe- 
riod. If this proves to be indeed the case, then theories that 
tie the Blazhko modulation to a side frequency of a pulsation 
frequency compo nent (kfo ± fm) m isinterpret the Blazhko 
phenomenon (e.g. iDziembowski fc M izcrski 2004) ■ 

• The global mean physical parameters of the star change 
about 1 — 2% during the Blazhko cycle. Most probably the 
primary changes occur in the mean stellar luminosity and 
the detected changes in mean stellar temperature, radius, 
and pulsation period are the consequences of the luminosity 
changes. 

• The phase modulation during the Blazhko cycle refiects 
simply the oscillations of the pulsation period, while the 
amplitude changes are due to the changes in the mean global 
physical parameters of the star (L, T, R) reflecting periodic 
alterations in the atmosphere structure of the star. 

Numerical models of the onset of radial nonlinear pul- 
sation display very similar dependence of mean radius and 
effective temperature on amplitude to the observed ones. It 
also suggests that, at least in part, the variation of mean 
quantities can be understood by the simple dynamics of 
radial pulsation. However, to unfold the whole picture, a 
presently unknown mechanism of the variation of the inter- 
nal structure of the star (e.g., turbulent/convective proper- 
ties) needs to be revealed. 
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Similar conclusions has been recen tly drawn abou t the 
nature of the Blazhko phenomenon bv lStothersI (|2006l l. He 
supposed that the underlying mechanism disturbing period- 
ically the amplitude and period of the purely radial mode 
pulsation is cyclic weakening and strengthening of turbu- 
lent convection in the stellar envelope. The suggested trig- 
gering mechanism, namely a turbulent convective dynamo, 
would, however, most probably result in much less regular 
modulation behaviour than observed in MW Lyr. Another 
difficulty with Stother's idea is how the multiperiodic na- 
ture of the modulation observed in some Blazhko stars e.g., 
CZ Lac (Sodor 200^ and Sodor, Jurcsik et al. in prepara- 
tion) and XZ Cyg ( LaCluvze et al.ll2004h can be explained. 
Consequently, while we think that though our results sup- 
port the idea that during the Blazhko cycle the pulsation 
remains purely radial, further theoretical efforts are needed 
in order to find the reason for the drastic changes observed 
in its main properties. 
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